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Abstract 

We consider a slow passage through a point of loss of stability. If the passage 
is sufficiently slow, the dynamics are controlled by additive random disturbances, 
even if they are extremely small. We derive expressions for the 'exit value' distri- 
bution when the parameter is explicitly a function of time and the dynamics are 
controlled by additive Gaussian noise. We derive a new expression for the small 
correction introduced if the noise is coloured (exponentially correlated). There is 
good agreement with results obtained from simulation of sample paths of the ap- 
propriate stochastic differential equations. Multiplicative noise does not produce 
noise-controlled dynamics in this fashion. 
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1. Introduction 

Bifurcations are normally observed by slowly changing a parameter ('sweeping') and 
waiting for a qualitative change in behaviour. The standard theory of bifurcations, 
however, does not have explicit time-dependence of parameters [1]. Here we study 
the effect of noise on the loss of stability in the simple equation 

y = ay- (i) 

This is a linearised dynamic bifurcation: a bifurcation in which the parameter g is 
explicitly a function of time. 

Normally it is said that loss of stability in (1) is signalled by exponential growth 
of |y| for (7 > 0. However if we continuously increase g from g <^ and wait for \y\ 
to exceed some fixed value, we find that this occurs well after ^ = [2]. In this 
paper we set 

9 = l^t (2) 
and let y take (fixed) initial value j/o at some ^ = < 0. Then the solution of (1), 

y = yoe^^^'-^°\ (3) 

reveals that \y\ < \yo\ until g = —go, regardless of the magnitude of go and how 
small the sweep rate fi is. However if n is sufficiently small the dynamics enter a 
noise-controlled regime, in which the value of g when \y\ reaches \yo\ is a random 
variable controlled by the size of additive random disturbances [3,4,5]. This random 
variable we call the exit value; it is more akin to a mean first passage time than to 
the usual definition of a bifurcation point [6]. 

This work is chiefly concerned with the probability distribution of the exit 
value in the "noise-controlled regime" [7] when (1) is perturbed by additive white 
and coloured noise. Because (1) is a non-autonomous linear equation, it is possible 
to calculate an exact solution for the probability distribution of y as a function of 
time and thus to calculate the probability that \y\ > \yo\, which we will denote 
by Pr(|j/| > \yo\). The probability distribution of the exit value is the derivative 
with respect to of Pr(|j/| > |yo|)- Additive noise reduces the average exit value by 
defining a level below which (y^) does not fall. 

The noise-controlled regime is in force if 

fi\lne\<^gi + 0{fi) (4) 



1 



where e is the (r.m.s.) noise level. This is the most likely regime for small fj,. Con- 
centrating on this regime makes possible a simpler analytical treatment than pre- 
viously reported [8, 6], including explicit derivation of the probability distribution 
of exit values. This is presented in section 2. We have chosen \y\ = \yo\ as the fixed 
exit level; this makes transparent the difference between the noise-controlled regime 
and the deterministic limit, but we emphasise that, as long as the noise-controlled 
regime remains in force, the dynamics are independent of the initial conditions. 

Our analytical results are in excellent agreement with results obtained from 
simulation of sample paths. These simulations we performed using an algorithm 
with weak second order convergence (the errors in the moments of the exit value dis- 
tribution are proportional to the second power of the timestep [9]). Our results are 
also consistent with the reported results of analogue and digital simulations [10, 6] 
performed in different parameter ranges. 

If the noise is coloured, the probability distribution of y is narrowed and the 
mean exit value is increased [8, 5]. The correction is, however, small except for very 
large values of the correlation time r [11]. In section 3 we derive a new formula 
for this correction which we successfully compare with numerical results in the 
(physically realistic) regime //r^ < 0{1). In section 4 we show that multiplicative 
noise has a qualitatively different (and in general much less dramatic) effect on 
the dynamics. In section 5 we discuss the situation where y has spatial degrees 
of freedom and the passage through the point of stability loss causes a growth of 
spatial order. 

Noise-controlled dynamics are useful in understanding the behaviour of sys- 
tems of ordinary differential equations used as low-order models in fluid mechan- 
ics [12, 13] and laser physics [14, 15]. In these cases the time-dependence of the 
parameter corresponding to g is in general more complicated; our system serves as 
a linearised model which, being a good approximation near g = 0, captures the 
essential dynamics. 
2. Additive white noise 

To study the effect of additive white noise we replace the ordinary differential equa- 
tion (1) by the stochastic differential equation (SDE) 

dyt = gytdt + edWt (5) 

where e is constant with < e -C 1 and Wt is the Wiener process [16]. The solution 
of this SDE with g — ij,t is 

Vt = e^^(*'-*o) (^y, + e j\-'^^(''-^'Uw,^ (6) 
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where to = go/fJ,. For each t, y is a, Gaussian random variable with mean equal 
to the solution (3) of the corresponding deterministic equation (1) and variance a 
function of time given by 



cr' = {yl) - {Vtf = e^e^'^ f e'^^^' ds. (7) 
If -go > 0{^) then, for g > 0{^), 



.-e^ei-y^. (8) 

The noise-controlled regime is in force when cr^ » (j/^) for g > 0{.^), so 
that all memory of the initial conditions {go, yo) is lost and noise controls the exit 
value distribution. For this we need 

/^Une| < i^o Qln^ -In^/o^ • (9) 

Since the right hand side is 0{1) in general, (9) holds even for extremely small 
values of e if /X is moderately small {fi < 0.1) [7]. 

In what follows, we assume that the noise level, e, is <^ \yo\ but large enough 
that (9) holds. For additive noise this means that, for g > 0{y/Jl), we can set 
{vt) — 0. The probability distribution of exit values is then given by 
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-Pr(|y|>|yo|) = ^^-- / e^dx 



yo 

2 I da f p r x"^ \ 

/ e^^dx + / —e^^dx I UUJ 



TT //cr^ dt 



2 gyo zlc 

e 2<^^ 

TT ji a 



using integration by parts. 

Let g be the value of g such that cr^(^) = j/q let v = —{g — g){g/ ji). Then 
cr ~ and the probability that the exit value of g lies between g and g + dg for 
additive white noise is RA{g)dg where, to lowest order in 



RA{g)^Jle-e^-'\ (11) 
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Figure 1. Distribution of exit values in the noise- controlled regime with white 

noise. Sample paths of (5) were simulated numerically numerous times, with 
1/0 = 1 and gQ = —1. The next value of g at which y = 1 was recorded 
each time and the resulting distributions of exit values are plotted along with 
the calculated distribution (11). The most probable value, given by (12), is 
indicated by a vertical line. (Without noise, the exit value is always 1; if 
fj,\ lne\ > 0.5 + (/x/4) /n(7r//i) the noise-controlled regime is not in force and the exit 
value distribution is narrow, Gaussian and centred on g = 1.) 



This distribution is compared with numerical results in Figure 1; it is the probability- 
distribution of ^ In where n is a Gaussian random variable with unit variance. 
The most probable value of v is 0, so the most probable exit value is g where 
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(12) 



We can obtain exact expressions for the moments of the exit value distribution 
in the noise-controlled regime by noting that f = + In 2) where v' — ln|n^. 
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The characteristic function, 0(s), of v' is known [17]: 



1 1 1 

= — / e-^^e^''e''''du = ^r(- + is). (13) 
V'^ J-oo V'^ 2 

The mean exit value, ((^e) is 

(^7e>=^ + ^(7 + ln2) (14) 

where 7 = 0.57721 . . . (Euler's constant) and and the variance of the distribution 
of exit values is 

(9?) - (9=)' = y fj. (15) 

It is instructive to consider how these results arise in terms of trajectories 
(sample paths) of (5). Firstly note that, once g > 0{^JJi)., most trajectories have 
\y\ » e, and so are no longer affected by the noise. Additive noise acts near ^ = 0, 
effectively providing a random initial condition for the subsequent evolution [18, 5]. 
The most probable exit value is determined by the time, 0(-\/ln(?/o/e)/A*), taken 
for the variance (j/^) to rise from a level proportional to at ^ = to j/q at 
g = g. The probability distribution of exit values has width C(/i) because most 
members of the ensemble of trajectories cross \y\ = \yo\ within a time 0{1) of each 
other, corresponding to an C(/i) spread in g. Now consider the exponential rise of 
neighbouring trajectories; the difference in their exit values is proportional to the 
log of the ratio of their starting points. The Gaussian distribution in y near ^ = 
is thus transformed into a distribution of exit values whose characteristic shape is 
that of I Inn^, where n is a Gaussian random variable. 

3. Additive coloured noise 

Noise is often added to an ordinary differential equation as a way of approximating 
the influence of neglected small effects while retaining simplicity of description. Of- 
ten "coloured noise" , i.e. noise with non-zero correlation time, is more appropriate 
than white noise. If y is subject to coloured noise, we replace (5) by 

dyt = ntytdt -|- er]tdt. (16) 

For rjt we use the most commonly studied coloured noise process [19], the Ornstein- 
Uhlenbeck process satisfying the SDE 

drjt = '^dt + dWt (17) 

T 
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which has 



{Vt)^0 and (ry^ry,) = i-e-l^-*!/- Vs, t. (18) 
The solution of (16) with y = 1 at t = to is 

= e^M(t'-*o) + e r e-'^(*'-*')r7(s)ds. (19) 

J to 

We restrict ourselves to the noise-controlled regime where, for t > 0{^), 

Jtn Jtn 



(20) 

where a = a//u/2(s — (/ur)~^) and cto = ^/Ji/2{to — (/ur)~^). Most of the support of 
the integrand is near s — 0. We therefore use the following approximation for the 
error function [20], 



/■a 
an 



e-^ ~ , (21) 

"0 — Q! + Ja'^ + ^ 

and expand the denominator in (21), about s = 0, in a form convenient for the 
integration over s: 

' 2 (22) 



where k = y^l + f A*'^^- Thus for t— r and |^o|~^ > (^(l/v^)' (j'l) approximately 
given by 

(.?>-;f.V.'e£J^. (23) 

This is smaller than the white noise value (8) by a factor which is a function of 
//r^. The solution of (16) also differs from that of (5) in multi-time correlations 
{ytVs) but this is not important for the purposes of calculating the distribution of 
exit values of g. We make the approximation that noise with correlation time r is 
equivalent to white noise at a slightly smaller level e' where 



e'~eJ^_eM-V8«^ (24) 
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The mean exit value in the presence of exponentially correlated noise with correla- 
tion time T is therefore given by: 



{g,) = y/2fi\\ne\-fififiT^) + + ln2) + ^(/x^) (25) 

where 

/(,.^)=ln(^)-£!-i.n(^), (26) 

7 = 0.57721 . . . (Eulcr's constant) and g is given by (12). This formula is compared 
with numerical results in Figure 2. The most probable exit value and variance of 
the exit value distribution can be calculated similarly. 

4. Multiplicative white noise 

If y is subject to noise which is proportional to y then we replace (1) by the SDE 



dy = gytdt + eytdWf (27) 

(The most natural way for multiplicative noise to arise is if g itself is subject to 
fluctuations.) The solution of (27) with y = yo at t — to is [16] 



y = yoe 



|M(t^-to)-|e'(*-*o)ge(m-mo) (28) 



where to = go/jJ-- The median value of y is j/oexp(^(5r^ ~ 9o ~ ^^{g ~ 9o)))', the 
factor exp(e(Wt — Wt^)) is always positive with mean exp(|e^(t — to)). As before, 
we first calculate the probability that y > yo as a, function of g. 

Pr{\y\ > \yo\) = Pr(e^(^*-^*o) > -) = / e ^dx (29) 

m V27r«T^J-lnm 

where m = exp(i^(t^ — tl) — |e^(t — to)) and cr^ — e^{t — to). When e ^ 1 the 
possibility that one trajectory crosses \y\ = \yo\ more than once can be ignored; 
the exit value distribution is then Gaussian, centred on — + with variance 
i2M)e'. 

Using our exit value definition, we find that multiplicative noise increases the 
median delay even beyond the deterministic value. Other definitions, based for 
example on the trajectory of [4, 8], may lead to the opposite conclusion because 



(^yt) and y are greater than the median value of y. Regardless of definition, 
the effect is small if e is small. For multiplicative noise, there is no analogue of 
the noise-controlled regime; we need 0{1) multiplicative noise to produce 0{1) 
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Figure 2. Distribution of exit values in the noise- controlled regime with coloured 
noise. In (a) and (b) the distribution obtained by simulation of sample paths 
solution of (16) is shown along with the calculated distribution (11) (smooth 
curve) calculated with e replaced by e' , according to (24), to account for the effect 
of non-zero correlation time in the noise. Graph (c) shows the calculated mean 
exit value as a function of r (25) with /x = 0.01 and e = 10~^^ (solid line) and 
some numerical results (crosses). 



effects. We have demonstrated this when the noise is simply proportional to y; it 
will remain true if the noise is any 0{1) function of y because the noise also becomes 
tiny as y does, and therefore does not dominate at any stage. In the presence of 
both additive and multiplicative noise, an exact solution can still be found for the 
probability distribution of y as a function of time; to derive the exit value statistics, 
it is useful to note that additive noise acts primarily near g = and multiplicative 
noise later [18, 11]. 
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5. Conclusion 

When the sweep through the point of stabihty loss is slow, the dynamics are con- 
trolled by very small random additive noise and memory of initial conditions is lost. 
For small values of the sweep rate (fj, < 0.1), noise produces 0{1) effects even if 
its magnitude is extremely small. In this noise-controlled regime the most probable 
value of g when y reaches an 0{1) predetermined level is ^ = \/2n\ lne| + C(^) 
and the probability distribution of exit values has standard deviation proportional 
to As long as the noise is additive and Gaussian, its correlation time affects the 
exit value distribution only slightly. Multiplicative noise does not produce noise- 
controlled dynamics in this way. 

A physical picture of this noise-controlled system is the dynamics of a strongly 
damped object allowed to slide in a slowly-changing potential and continuously 
subject to random influences. The variable y is the distance from the extremum 
of the potential. For t < 0, the extremum is the centre of a well; for t > the 
extremum is the crest of a hill. Once t > 0, the object lingers near the crest for a 
time controlled by the noise level before sliding away. 
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